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ABSTRACT 

In a previous paper we demonstrated that non-radial hydrodynamic oscillations of a thermally- 
supported (Bonnor-Ebert) sphere embedded in a low-density, high-temperature medium persist for 
many periods. The predicted column density variations and molecular spectral line profiles are similar 
to those observed in the Bok globule B68 suggesting that the motions in some starless cores may be os- 
cillating perturbations on a thermally supported equilibrium structure. Such oscillations can produce 
molecular line maps which mimic rotation, collapse or expansion, and thus could make determining 
the dynamical state from such observations alone difficult. 

However, while B68 is embedded in a very hot, low-density medium, many starless cores are not, 
having interior/exterior density contrasts closer to unity. In this paper we investigate the oscillation 
damping rate as a function of the exterior density. For concreteness we use the same interior model 
employed in Broderick et al. (2007), with varying models for the exterior gas. We also develop a 
simple analytical formalism, based upon the linear perturbation analysis of the oscillations, which 
predicts the contribution to the damping rates due to the excitation of sound waves in the external 
medium. We find that the damping rate of oscillations on globules in dense molecular environments 
is always many periods, corresponding to hundreds of thousands of years, and persisting over the 
inferred lifetimes of the globules. 
Subject headings: stars: formation, hydrodynamics, ISM: globules, ISM: clouds 



1. INTRODUCTION 

The small dark molecular clouds known as starless 
cores are are significant in the interstellar medium 
as the potential bi rthplaces of stars (review by 
iBergin fc Tafallall2007[ ). As their name implies, the star- 
less cores do not yet contain stars, but their proper- 
ties are nearly the same a s similar small clouds t hat do 
(iMvers. Linke. fc Benso"^ [l983: Mvers fc Bensoni 119831 
iBenson fc Mverg Il989f) . Furthermore, there is a com- 
pelling similarity in the mass function of the starless cores 
(jLada et al.ll2008l ) and the initial mass function (IMF) of 
stars. 

Observations of Bok globules (Bok 1948) 

and starless cores (|Ward -Thompson ct al. 

1991 iTafalla et al. I [l998l: llee fc Mvera ,199a: 
Lee. Mvers. fc Tafallal 120011: lAlves. Lada. fc Ladal 
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Bacman: 
200l[ 



Ward-Thompson, Motte, fc Andre 1999; 

i_et al. 200^ Shirlev et al. 2000; Evan s et al. , 
Shirelv ct all [200l lYoung et al. I l200l 
Tafalla e~al. „2004i : lKeto et al.|[2004D suggest that many 
of these small (M < 10 M©) clouds are well described as 
quasi-equilibrium structures supported mostly by ther- 
mal pres s ure, a pproximately Bonnor-Ebcrt (BE) spheres 
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Aguti et al.l l2007f) 



show complex profiles that further suggest velocity 



^ Current address: Canadian Institute for Theoretical Astro- 
physics, 60 St. George St., Toronto, ON M5S 3H8, Canada 



and density perturbations within these cores. We have 
previously shown that in at least one case, B68, these 
profiles can be produc ed by the non-rad i al oscillations of 
an is othermal sphere (jKeto et al.l 120061 : iBroderick et al.l 
[20071) . 

To match the spectral line profiles, we have previously 
found that the oscillations have to be large enough (2 5 
%) that the amplitudes are non- linear (iKeto et al.ll2006f ). 
If such large amplitude oscillations are to be a viable 
explanation for the observed complex velocity patterns 
then their decay rate must be no faster than the sound- 
crossing or free-fall times of the globules, 10^~^ yrs. Nu- 
merical simulations in which the external medium is sub- 
stantially less dense than the gas inside the core have 
found that mode damping is s ufficiently slow that m odes 
will persist for many periods ([Broderick et al.ll2007f ). 

However, it is not clear that this description is ap- 
propriate for many starless cores. B68 is one of only a 
handful of cores surrounded by hot, rarefied gas in the 
Pipe Nebula, where the vast majority of cores, like those 
in the Taurus fc Perseus molecular clo uds, appear to be 
embedded in cold , dense molecular gas ([Lada et al.ll2008l : 
lGoldsmith1l2008f) . Embedded oscillating clouds will gen- 
erally act as sonic transducers, exciting sound waves in 
the external medium and thereby loosing energy. A sim- 
ple one-dimensional analysis, discussed in detail in the 
Appendices, suggests that when the density contrast be- 
tween the core interior and its bounding medium is close 
to unity this process can be very efficient, potentially 
limiting the lifetime of oscillations. 

In this paper we investigate the damping rate of oscil- 
lations of isothermal spheres as a function of the density 
contrast with the external bounding medium. We do 
this using numerical simulations of a large amplitude os- 
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Fig. 1. — The initial, unperturbed density profile of the embed- 
ded isothermal sphere. The density contrast for each is roughly 50 
(black), 25 (blue), 10 (green) and 4 (red). This color scheme will 
be consistent throughout the other figures. 



TABLE 1 
External Gas Properties and Damping Rates 



Model Text/T pext/p.urf" TA [(GM/B?)- 
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green 
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red 
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3.44 
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- 1.24 X 10^ yr 





dilation superimposed upon an isothermal sphere. To 
evaluate our results we also derive analytical and semi- 
analytical estimates of the damping rate associated with 
the excitation of sound waves in the exterior medium. A 
brief discussion of the numerical simulation is discussed 
in section [21 the presentation and discussion of the nu- 
merical results are in section [31 and conclusions can be 
found in section[31 The details of the analytical and semi- 
analytical computation of the damping rates of small am- 
phtude oscillations are relegated to the appendices. 

We find that a reduced density contrast between the 
core and the exterior bounding medium does result in 
increased dissipation of oscillations as momentum and 
energy are transferred out of the core through the bound- 
ary. However, if the oscillations are non-radial the dis- 
sipation through the boundary is always less than the 
dissipation rate due to mode-mode coupling. 

2. NUMERICAL HYDRODYNAMIC MODEL 

Because the thermal gas heating and cooling time, 
via coUisional coupl i ng t o dust at high densities 
(jBurke fc HoUenbach I Il983f) and molecul ar line radia- 
tion at low densities (Go ldsmith I |2001[ ). is short in 
comparison to the typical oscillation period (on the 
order of 10^ yr), the oscillations of dark-starless cores 



are well mode l ed wi th an isothermal equation of state 
(jKeto fc Fieldl l2005| ). In our simulations, we use a 
barotropic equation of state with an adiabatic index of 
unity (as opposed to 5/3, for example), in which case 
the isothermal evolution is also adiabatic. This allows us 
to replace the energy equation with a considerably sim- 
pler adiabatic constraint. Since the unbounded isother- 
mal gas sphere is unstable, we truncate the solution at 
a given center-to-edge density contrast {p{R)/pc ~ 12, 
not to be confused with the interior/exterior density con- 
trast), producing a stable Bonnor-Ebert sphere solution 
with radius (R) and mass determined by the central den- 
sity (pc) and temperature (T). This is done by inducing 
a phase change in the equation of state, i.e.. 



P = 



p(R) 



P^ iip>p{R), 
PiR)^ ifp(i?)>P>/3p(i?)5^, 



-J "--^ cxt 
filTlp 



otherwise , 



(1) 

where the intermediate state, with /3 = 0.8 and e = 
log(/?)/log(/3r/Tcxt) < 1 chosen to make P{p) contin- 
uous, is introduced to avoid numerical artifacts at the 
surface, and Text ^ T. The specific values of Text that 
we chose, their corresponding density contrasts and de- 
cay timescales, are presented in Table [1] The associated 
radial density profiles are shown in Figure [H Note that 
there is a transition region which distributes the sud- 
den drop in density over many grid zones, though the 
pressure is continuous (necessarily) and nearly constant 
outside the surface of the Bonnor-Ebert sphere. 

We employ the same three-dimensio nal, self- 
xavit ating hydrodynamics code as in Broderick et al.l 
20071 ) to follow the long term mode evolution. Details 
regarding the numerical algorithm and specific valida- 
tion for oscillating gas sp heres (in that case a whit e 
dwarf) can be found in IBroderick fc Rathore I (120061 ). 
and thus will only be briefly summarized here. 

The code is a second-order accurate (in space and time) 
Eulerian finite-difference code, and has been demon- 
strated to have a low diffusivity. Because the equation of 
state is barotropic, we use the gradient of the enthalpy 
instead of the pressure in the Euler equation since this 
prov ides better stability in the unperturbed configura- 
tion ([Broderick fc Rathore |[2006r ). The Poisson equation 
is solved via spectral methods, with boundary conditions 
set by a multipole expansion of the matter on the compu- 
tational domain. As described in Broderick et al. (2007), 
the code was validated for this particular problem finding 
convergence by a resolution of 128'^ grid zones. 

The initial perturbed states are computed in the linear 
approximation using the standard form alism of adiabatic 
(which in this case is also iso thermal, iKeto et al1l2006[ ) 
stellar oscillations ([Coxlll980l ). The initial conditions of 
the perturbed gas sphere are then given by the sum of 
the equilibrium state and the linear perturbation. We 
employ the same definition for the dimensionless mode 
amplitude as given in the Appendix of Kcto ct al. (2006). 

Throughout the evolution, the amplitude of a given 
oscillation mode, denoted by its radial (n) and angular 
(/ fc m) quantum numbers, may be estimated by the 
explicit integral: 



-^nlm — 



.L/d^-pv.cL„, 



(2) 



where v is the gas velocity, ^nim i^ the mode displace- 
ment eigenfunction, uj„im is the mode frequency and the 
dagger denotes Hermitian conjugation.^ 

3. RESULTS & DISCUSSION 

In principle, a number of hydrodynamic mechanisms 
exist by which the oscillations can be damped. In the 
absence of a coupling to an external medium, these are 
dominated by nonlinear mode-mode coupling, in which 
large scale motions e xcite smaller scale perturbations 
(jBroderick et al.ll2007f l. However, when embedded in a 
dense external medium it is possible for the mode to 
damp by exciting motions in the exterior gas. That is, 
it is possible for the oscillating sphere to act as a trans- 
ducer, generating sound waves in the external gas which 
then propagate outward resulting in a net radial energy 
flux and damping the oscillations. 

Generally, this will be a strong function of the density 
contrast between the interior of the core and the exte- 
rior medium. If the external density has a low density, 
and thus little inertia, outwardly moving sound waves 
will contain little energy density despite their large am- 
plitude and thus inefhciently damp the mode (as is the 
case for B68). Similarly, if the exterior medium has a 
very high density the excited outgoing waves will have 
small amplitudes, which despite the large gas inertia will 
again carry away only a small amount of energy. Con- 
versely, when the density contrast is nearly unity, there 
will be an efficient coupling between the interior and ex- 
terior waves, resulting in a rapid flow of energy out of 
the cloud. 

This may be made explicit via the standard three- wave 
analysis, in which propagating wave solutions for the in- 
cident, reflected and transmitted waves on either side of 
the density continuity are inserted into the continuity 
and Euler equations, which can subsequently be solved 
for the reflected and transmitted energy flux. In this ide- 
alization the ratio of the reflected to transmitted flux is 
4C/(1 + C)^ where 



C = PeCs,e/ PiCsA = \/ Pel Pi 



(3) 



is the ratio of the sonic impedances on either side of the 
discontinuity ( Landau fc Lifshitz Ill987t l. As anticipated, 
when the density of the external gas is comparable to the 
surface density of the isothermal sphere we may expect 
efficient transmission of sound waves from the interior to 
the exterior (i.e., conversion of the oscillation into travel- 
ing waves in the exterior), and therefore efficient damping 
of the pulsations. 

This conclusion is also supported by a ID analysis of 
the decay of a standing wave confined to a high density 
region (section IA.2p . for which the damping rate, 7 is 
given by 



2R 



In 



li-CI 
1 + C 






(4) 



where Cs^i is the sound speed in the interior and i? is 
the radius of the high-density region. Thus, as expected 

■^ For stellar oscillations the mode amplitudes may be determined 
from the density perturbation alone using the velocity potential. 
However, when there is a non-vanishing surface density, as is the 
case for the Bonnor-Ebert sphere, this is not possible and the ve- 
locity field must be used. 



lower external densities (corresponding to smaller C) 
have lower damping rates and larger damping timescales. 
Note that when C is of order unity the wave decays on 
the sound crossing time of the dense region, as would 
be expected if, e.g., we were to decompose the standing 
wave into traveling waves. 

More directly applicable to the damping of oscillating 
gas spheres is the case of the decay of a standing multi- 
polar sound wave in a uniform density sphere, discussed 
in section [A.3l This is distinct from the ID case in an im- 
portant way: the oscillation wave vectors are no longer 
only normal to the boundary. As a consequence, the 
damping rate is a strong function of not only C but also 
the multipole structure of the underlying oscillation. In 
particular, the damping rate for the Z**^ multipole mode 
is pr oportional to (^^(+3 ^^j. gj-^all ^ (see, eqs. IA23I & 
IB17p . which is a very strong function of ^ for even low 
l\ This limiting behavior is also found in semi- analytic 
calculations (section |B|) in which the linear oscillation 
of an isothermal sphere is properly matched to outgoing 
sound waves in the exterior. The damping timescales 
(inverse of the damping rates) from both the uniform- 
sphere approximation and the full linear mode analysis 
are shown in Figure [5] for the first few multipoles. All of 
this implies that the damping rates of oscillations of cores 
embedded in cold molecular regions, where the density 
contrast between the core and the cloud is near unity, 
will be much more rapid than those in cores isolated in 
hot, lower-density regions. 

However, this is not borne out by numerical simula- 
tions of high amplitude oscillations. This can readily be 
seen by the damping timescales measured via the nu- 
merical simulations. Figure [3] shows the evolution of the 
n = 1 , / = 2, m = 2 mode discussed in iBroderick et all 
(|2007f ). In all cases the initial amplitude was 0.25, and 
the evolution of the mode well flt by a decaying expo- 
nential. Oscillations on isothermal cores embedded in 
higher density regions did indeed damp more rapidly, 
as is apparent in the figure. Explicit values of the de- 
cay timescale a re presented in Table [1] though as in 
IBroderick et al.l (|20 07) these should be seen as lower 
limits only due to coupling in the artificial atmosphere. 
These are plotted as a function of the surface density con- 
trast in Figure [21 shown by the colored filled circles. For 
reference, we show th e ap proximate surface density con- 
trasts of Barnard 68 (lAJv cs. Lada, & Lada 2001), L1544 
(jWard-Thompson. Motte fc Andrelll9 99) and typical for 
cores in the Pipe Nebula (jLada et al.ll2008. ) . While there 
is a clear power- law dependence of the damping timescale 
upon C, with TA oc {pe/ Pi)~^'^'' this dependence is consid- 
erably weaker than that predicted by the linear analysis 

[ta oc (pe/p«)"^-^). 

This disparity may be understood in terms of the rel- 
ative importance of the excitation of sound waves in the 
exterior medium to the damping of the oscillations. In 
particular, it is notable that the numerically measured 
damping timescales are always less, and for most surface 
density contrasts considerably so, than those implied by 
the linearized analysis of external sound wave excitation. 
This is true even with low amplitude oscillations (e.g., 
A ^ 10^^, for which the damping timescale is roughly a 
factor of 3 larger than for A = 0.25 ) and thus is likely 
due to mode coupling in the transition region immedi- 
ately outside the cloud (seen in Fig. [T] immediately out- 
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Fig. 2. — The damping timescales in units of the dynamical time 
of the cloud, {GM/R^)~^'^, as a function of the ratio of the exte- 
rior and surface densities. For reference the corresponding damping 
times in Myr is shown on the right-hand axis for a cloud similar to 
Barnard 68. The solid lines show the damping timescale predicted 
by a linear mode analysis of the isothermal sphere for monopole 
(Oil), dipole (111) and quadrupole (122, thick) modes (see the ap- 
pendices). The asymptotic behavior of the linear mode analysis is 
shown by the short dashed lines for each, which may be obtained 
by treating the oscillations of the isothermal sphere as sound waves 
in a uniform density sphere (see the appendices). The solid circles 
show the damping timescales measured via numerical simulations, 
and are colored to match the curves in Figures [T] & |3] Finally, we 
show a power-law fit to the measured damping timescales, which 
have a power-law index of —0.37, and is considerably flatter than we 
might have expected (though these represent lower limits upon the 
decay timescale). For comparison the density contrasts observed 
in some well known examples are also shown, including those for 
Barnard 68 (10"^), L1544 (10"^) and the cores in the Pipe Nebula 
(0.1-0.7). 

side r/ R = 1 as the region of rapidly decreasing density). 
Within this region the densities are all within a single 
order of magnitude despite their very different asymp- 
totic densities at infinity (which differ by many orders 
of magnitude) . As a consequence, the damping rates in- 
ferred from the isolated isothermal spheres are only a 
weak function of the properties of the external medium. 
Thus, even for surface density contrasts on the order of 
1/4, the decay timescale for the quadrupolar oscillation is 
still larger than 4 x 10^ yr, corresponding to many oscilla- 
tions and comparable to the inferred lifetimes of globules. 

4. CONCLUSIONS 

Despite the expectations of the linear mode analy- 
sis of the embedded isothermal sphere, the damping 
timescale of large-amplitude oscillations of embedded 
Bonnor-Ebert spheres is only a weak function of the 
density of the external medium. This is a result of 
the dominance of nonlinear mode-mode coupling in the 
damping of large oscillations. Even in cold molecular 
environments, the quadrupolar oscillation discussed in 



iKeto et all (|2006D and lBroderick et all (|2007f ) has a hfe- 
time of at least 0.4 Myr, and i s thus comparabl e to the 
inferred lifetimes of globules (|Lada et al.l l2008f ) . This 
suggests that globules supporting large-amplitude oscil- 
lations may be common even in these environments. 

The presence of large- amplitude oscillations on starless 
cores has a number of consequences for the physical in- 
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t (Myr) 

Fig. 3. — The mode evolution for various external gas models. 
The density contrast for each is roughly 50 (black), 25 (blue), 10 
(green) and 4 (red), and the lines are colored to match those in Figs. 
[l] & [2] The decay time for the red curve is 0.43 Myr, compared 
to 1.05 Myr for the black curve. Quadratic fits for each are shown 
by the dashed lines. Matching these up to specific models of the 
mode decay is complicated by (i) the non-linear mode dynamics, 
and (ii) understanding the energy losses due to traveling waves in 
the external medium (see the appendices). 



terpretation of observations of starless cores (jKeto et al.l 
120061: iBroderick et al.ir2007[ ). Necessarily they would im- 
ply that starless cores are stable objects, existing for 
many sound-crossing times. In addition, the signatures 
of collapse, expansion and rotation in observations of self- 
absorbed molecular lines (e.g., CS) are degenerate with 
molecular line profiles produced by oscillating globules 
(see, e.g., figure 6 of Broderick et al. 2007). Conse- 
quently, it may be difficult to determine the true dynam- 
ical nature of motions in a given globule or dense core 
from self-absorbed, molecular-line profiles alone. This 
suggests that studies of such profiles in dense cores may 
be of only statistical value in determining the general 
status of motions in dense core populations. 



We would like to thank Mark Birkinshaw for bringing 
the problem of damping starless core pulsations in dense 
media to our attention. 



APPENDIX 
DAMPING OF AN OSCILLATING CAVITY 



We will first review some basic facts about sound waves in uniform media. Beginning in subsection lA.2l we will discuss 
the evolution of a one-dimensional standing sound wave as a result of the excitation of waves in the external medium. 
In subsection I A. 31 we discuss the application to oscillating uniform density spheres, and treat the full linearized mode 
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analysis of isothermal spheres in section |B] 

Equations of Motion 
The governing equations are the linearized continuity and the Euler equations; 

p' + V • podw = , 

Sv+^Vp'^0, (Al) 

Po 

where po and Cs are the unperturbed density and sound speed, respectively, p' is the Eulerian perturbation in the 
density and dv is the Lagrangian perturbation in the velocity (which is identical to the Eulerian perturbation since the 
initial velocity field is assumed to vanish) . These may be combined in the normal way to produce the wave equation 

p'-clV^p' = 0, (A2) 

where we have assumed that the background density is uniform. Solutions will generally obey the dispersion relation 
k^ = <?gUp' ^ and may then be inserted into the Euler equation to obtain the associated velocity perturbation. 

Decay of a One- Dimensional Sound Wave 

We will consider the decay of a sound wave traveling inside of a uniform high density region (called the interior and 
denoted by sub-script i's) surrounded by a uniform low density region (called the exterior and denoted by sub-script 
e's). The interior wave oscillation is given by 

p'. = Aipi sm{kix)e~'''^* and Svf = -iAiCs,i cos(fcia;)e"'"* . (A3) 

The exterior waves are outgoing, and we will focus upon the boundary condition at the +R boundary. At this point 
the exterior sound wave is given by 

p'^ = A^p^e'''"''"'^' and 6v^ = AeC,,ee'''='=~'"* . (A4) 

At the interface we require (i) pressure equilibrium and (ii) continuity of the displacement and hence velocity. Since we 
have assumed the pi and pe are constant, the first condition is simply p'^c^ ^ = p'^c^ g. The second is trivially Svf = Sv^. 
Thus 

A,p,cl, sm{hR)e-"^' = AePee'''''^-"^* and - ^Acs,^ cos(fc.i?)e-'"* = AeC.,ee'''=-"-*"* , (A5) 

give two complex equations from which we may determine Ae and lo as functions of Ai and the background fluid 
quantities (note that ki = uj/cs,i and ke = uj/cs,e by the dispersion relations in each region). In particular, note that 
by taking the ratio of the two equations we find 

tan(fc,i?) = -< , (A6) 



where C, = PeCg.e/ PiCs.i — \/ Pel Pi is the ratio of the exterior and interior sonic impedances. Generally this may be 
solved to find that kiR = riTr -f (i/2) In [|1 — C|/(l + C)]; however we will solve this explicitly in the C ^ 1 limit to 
illustrate how we will do this for spherical geometries later. 

Let us begin by assuming that the damping rate is small. Specifically, let us assume that loq = 5R(aj) is much larger 
than 7 = — 3(cj). Then we may Taylor expand the left-hand side of eq. (jA6p around 7 = 0: 

<, (A7) 

and thus, equating real and imaginary parts 

(A8) 

and 

(A9) 

We may in principle then insert this into eq. (jA6p to obtain the relationship between Ai and A^. However, we will be 
concerned with only the damping rate here. 

Damping Timescale of an Oscillating Sphere 

In the previous section we discussed the simple problem of the damping of a one-dimensional wave due to the 
excitation of outgoing exterior sound waves. In this section we will address the more relevant problem of the damping 
timescale of multipolar oscillations on a sphere. While generally we would use the linearized analysis of the oscillating 
isothermal sphere, here we will limit ourselves to the spherical analog of the previous section: a sound wave in a 
uniform density sphere surrounded by a uniform density exterior. 
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Before we can discuss the excitation of sound waves outside of an oscillating sphere we must first determine the 
explicit form of these waves. We do this by separating the radial and angular dependencies using spherical harmonics 
(primarily because these were used in determining the mode spectrum of the isothermal sphere). That is, we let 
p' — e^'^'^*TZ{r)Yim,{i'), and insert this into equation (|A2p . producing 

l^rn^^-^^^^n + k'n^o, (aio) 

where k = lo/cs- The general solutions of this equation are simply the spherical Bessel functions: 

TZ{r) = aji {kr) + bni {kr) . (All) 

However, we must still separate the inward and outward traveling waves. This naturally occurs if we consider spherical 
Bessel functions of the 3rd kind, which are necessarily simply linear combinations of spherical Bessel functions of the 
1st and 2nd kind (j/ and nj, respectively): 



K' i^) =jiiz)±ini{z) oc -jq^Y 



±22 



(A12) 



,(1)^ 



,(2), 



Outwardly directed waves are given by h^ (z) and inwardly directed waves are given by h\ (z). In terms of this the 
density perturbation of the exterior sound waves associated with the l,m multipole is 



p'^ = A,p,e-''^'h\'\ker)Yi^ir) 



(A13) 



where A^ is a dimensionless wave amplitude. For I = this results in the standard spherical wave solutions, oc e^*'^''"/r. 
The velocity perturbation may be determined via the linearized Euler equation: 



Sv: 



ujpe Or 



,(1) 



Z/leCo_e6 



--*^''' ■(Ze)>l™(f), 



dZe 



(AM) 



where here and henceforth we have defined Ze = k^r for convenience. In contrast, the interior wave must be regular at 
the origin, and is given by 



p[ = A,p,e-''^'ji{hr)Yi„,{v) and 5< = -^A;c,,,e-*'^*■^(z,)r^™(f) , 

dz^ 

where Zi = kir = Ze/C- 
At the interface we again require pressure equilibrium and continuity in the radial velocity: 



2 ^-i^tAi), 



A,p,ci^^e-'^'ji{hr)Yi^{v) = A^p^ci^^e-'^^'hY' {k,r)Yi^{v) 



and 



UlJ 
lAiCs iS — — (Zj)y;„i(r) . = —ICse^ 

ozi 



-iut I 



dh 



(1) 



dZe 



-{Ze)Yim{v) 



(A15) 



(A16) 



(A17) 



As before, together with lj — ki/ Cs,i — ke/cs.e and z^ — (zi, these are sufficient to determine A^ and lu. In particular, 
these imply 

dlnji 



dzi 



-(zi 



{Ze) 



C dZe 



(A18) 



The right-hand side may be simplified in the Zg <^ 1 limit. Since Zi is typically of order unity, this implies that C <C 1. 
In this regime, to leading order in small Ze, 



Idlnh 



(1) 



C dze 



<Ze) 0^ 



l + l 



+ i- 



z [(2/-1)!!]' 



(A19) 



where {21 - 1)!! = (2^ - 1) ■ (2^ - 3) • ■ • 5 ■ 3 • 1. 
Again let us begin with the ansatz that the damping rate is small for small C. In which case 



dlnji 



dz,. 



d'^lnji 



.jR 



l + l 



2l~l.,2l 



C~'z. 



Cz, [(21 - 1)!!] 



(A20) 



and thus 



dlnji 



jl~l{Zi) _ /+ 1 

3i{zi) Zi 



l + l 

jl-lizi) l + l J Cs,i „ 



(A21) 



Embedded Oscillating Starless Cores 



where Zni is the n root oi ji{z) 



With 



dzf 



Jl-l{Zi) 



3i{zt) 






(where we used the fact that C, <^\) we find 
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C 



21+3 ■y'21+2 



R {I + 1)^(21 - ly.iy 



(A22) 



(A23) 



We note that this is quite different than the one-dimensional case. In particular, the decay rate decreases rapidly 
with decreasing (^, and is a strong function of I, with higher multipoles decaying considerably more rapidly than lower 
multipoles. This is a consequence of our assumption that not just C <C 1, but (^Zi <C 1 (justified in our case), which 
implies that Ae ^ R. That is, the exterior propagating sound wave necessarily knows that the geometry is converging, 
with higher multipoles converging more rapidly. 

We should also emphasize that the scaling of the damping rate, 7, with C is dependent primarily upon the structure 
of the traveling sound waves in the exterior. Rather, the properties of the interior perturbation are encoded in 
the definition of the Zni, and thus the normalization. Therefore, the damping rates associated with the linearized 
perturbations of the pressure-supported isothermal sphere (Bonnor-Ebert sphere) should be quite similar, differing 
only in the particular values of Zni — 27ri?/Aj. Indeed, in the next section we find this to be the case. 

DAMPING OF OSCILLATING ISOTHERMAL SPHERES IN THE LINEAR REGIME 

The rather surprisingly strong scaling of the damping timescale with the density contrast motivates a more careful 
analysis of the damping rates of a self-gravitating, pressure supported isot hermal sphere. The mode a nalysis of 
an isothermal sphere has been treated in considerable detail elsewhere (e.g., ICoxJll980t iKeto et aD 120061 ) and thus 
we summ arize it only briefly here. The perturbed quantities are most conveniently described by the Dziembowski 
variables ()Dziembowski|[l97ll ) 

V "^ Pogr gr g or J 



which may be separated into radial and angular parts, ?7i(r) 
hydrodynamic equations are given by 



77,(r)e-*"*y^ 



;(f). In terms of these, the linearized 



dr ^ 



1(1 + 1) 



dr]2 
dr 
dri3 

dr 

dri4 

dr 



where 



U 



dlnr^g 
dlnr 



V 



3)^1 , ^ ^,^ 
= a^Cr]i + (1 - C/)772 
= (1 - [/)?73 + V4 
^ UVrj2 + [l{l + 1) - 

dlnPo 



V 



m + Vr]^ 



UV] 773 - C/ry4 
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C 



GM 



& 



3, ,2 



GM 



(B2) 

(B3) 
(B4) 
(B5) 

(B6) 



are functions of the background equilibrium configuration. These are solved subject to a normalization condition 
771 (i?) = 1 and the inner boundary conditions: 



f^ C'7?i „ = /ryz L ^n 



and 



(B7) 



arising from regularity at the origin. A third boundary condition is a result from the continuity of the gravitational 
potential at the sphere's surface, 

ml^R = -{i + ^)ml^n- (B8) 

If the exterior density vanishes, as was assumed in lKeto et al.l ()2006l ). the final boundary condition is given by requiring 
that the Lagrangian pressure perturbation, 

5P = P' + 8r- VPo = pogr [m - V3 - Vi) ^ ■ (B9) 

However, we now properly match this onto the outgoing multipolar wave solutions described in section [A. 3 1 Namely, 



6P ^ p'.cl^ ^ A,p,h\^\z,) , 



subject to 



5r = rjiR 



Uj'^Pe 



dr 



A, 



OJ 



dzp^ 



(1) 



(BIO) 



(Bll) 



Therefore, 



m -V3 -Vl\r=R = Vl 



Broderick et al. 



r=R 



dhW/dze pog 
As before we assume that the damping rate is small, set to = luq — ij and employ the boundary condition 



m -m -m\r=R ^'^ 

The damping rate is then estimated by 



m 



dhW/dze Pog 



r—R,uj—LJo 



7 = — Q^ 



/i(i) 



Vi 



PeCs,e^ 



dhW/dze Pag 



dim - m) 



r—R,Lj—LUQ • 



du 



(B12) 



(B13) 



(B14) 



r=R 



where we used the fact that drji/dcu = by the normalization condition. In practice 9(772 — r]^)/duj is evaluated by 
solving the eigenvalue problem for the ry^ and a^ using 



m -m -'ni\r=R = "^ 



/id 



Vi 



PeCs,e^ 



dhW/dze Pog 



dim - m) 



r—R.^uj—ujo . 



duj 



r=R 



and setting 



dim ~ m) dim ~ m)/dc 



did 



dbj/da 



(B15) 



(B16) 



This amounts to finding dim — m)/d'jJ holding the other boundary conditions constant. 

The results of this procedure are explicitly shown in Figure [51 Of particular note is that the intuition obtained from 
the analysis of sound waves in a uniform density sphere is borne out in the low pe/ poiR) limit, for which 



7 ex 



(2/+3)/2 
/^e \ _ /-2/+3 



PoiR) 



(B17) 
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